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Some commonly used approximations for complete fusion and breakup trans- 
mission coefficients in collisions of weakly bound projectiles at near barrier 
energies are assessed. We show that they strongly depend on the adopted 
classical trajectory and can be significantly improved with proper treatment 
of the incident and emergent currents in the WKB approximation. 
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I. INTRODUCTION 

In a fusion reaction the projectile and target nuclei form an excited compound system, 
which decays by emission of particles or gamma rays. When there are channels strongly 
coupled to the elastic channel, the reaction may be described using any of the available 
coupled channels codes. The effects of such couplings have been extensively discussed in 
the literature flj . In general terms, at energies below that of the Coulomb barrier, these 
couplings tend to reduce the effective fusion barrier, which substantially increases the fusion 
cross section. At high energies the incident flux in the elastic channel is partially diverted into 
inelastic and transfer channels. This tends to decrease the contribution of the elastic channel 
to the fusion cross section, but this reduction is partially compensated by the contribution 
from the other channels. 

The recent availability of radioactive beams has made possible to study reactions involv- 
ing unstable nuclei. Such reactions are important in processes of astrophysical interest, as 
well as in the search for superheavy elements. The main new ingredient in reactions induced 
by unstable projectiles is the strong influence of the breakup channel. In the case of not 
too unstable projectiles, the effect of this channel in the fusion cross section at low energies 
is, as in the case of stable beams, to enhance it. At high energies, however, the situation is 
qualitatively different from the case where only stable nuclei are involved. The contribution 
from the breakup channel to the fusion reaction is strongly influenced by the low probability 
that all fragments are captured. Thus, in this case, the fusion cross section is partitioned 
into a complete and one or more incomplete fusion contributions [0. 

The introduction of the breakup channel into a coupled channels calculation is by no 
means trivial. The difficulty lies in the fact that this channel lies in the continuum, and 
involves a, at least, three body system. This problem has been addressed by several authors, 
using different approaches. Several recent experiments involving fusion of neutron rich 6 He 
and proton rich 17 F with heavy targets have been performed with the purpose of exploring 
these theoretical proposals |3|-|5| . In Refs. 0^-0 , the coupled-channel problem is simplified 
by the introduction of the polarization potentials arising from the coupling with the breakup 
channel [pTG| , pm[j . In Refs. |T^-|n|, the coupled-channel problem is solved directly within 
different approximation, ranging from the schematic model of Dasso and Vituri [12| to the 



huge calculation of Hagino et al. |T3[], performed through continuum discretization. 

The polarization potential approach of Refs. B,0-|9j has the advantage of leading to sim- 



ple expressions, which can easily be used in data analysis |15 |. However, it employs several 
approximations which were not thoroughly tested. These approximations can be grouped 
in two cathegories. In the first are those used in the derivation of the polarization poten- 
tials. In the second are the semiclassical approximations for fusion and breakup coefficients, 
used in calculations of the cross sections. These coefficients are written in terms of barrier 
penetration factors and survival probabilities, which are evaluated within the WKB approx- 
imation. The aim of the present work is to ascertain the quality of the approximations for 
the transmission coefficients. Approximations in the derivation of the polarization potential 
will be the object of a latter study. For our purposes, we consider a case where a complete 
quantum mechanical calculation is feasible and compare exact and approximated cross sec- 
tions. We study the n Li + 12 C collision, using typical optical and polarization potentials. 
For simplicity, our polarization potential has no angular momentum or energy dependence 



and the range is given by the n Li breakup threshold energy. The strength is consistent with 
that found in Ref. [JHJ for the most relevant partial waves in near-barrier fusion. 

The plan for this paper is as follows: in section II we briefly revise the coupled channels 
formalism and the concept of a polarization potential. In section III we discuss different 
approximation for the transmission coefficients and investigate their consequences on the 
fusion and breakup cross sections. Finally, in section IV, we present the conclusions of this 
work. 



II. COUPLED CHANNEL EQUATIONS AND POLARIZATION POTENTIALS 

In a standard coupled channels calculation, the system is described through the distance 
between centers of projectile and target, r, and a set of intrinsic coordinates, £, that describe 
the internal degrees of freedom of one of the nuclei, e.g. the target. These coordinates are 
associated to an intrinsic Hamiltonian h and its eigenfunction set, 

hMO = e«0a(O , (1) 

where 

f fc(t)MtW = 8** • (2) 

The system Hamiltonian may then be written as 

H = T + U opt + h + v(r,Z). (3) 

Above, T is the kinetic energy of the relative motion, U opt is the optical potential and v(r, £) 
is the interaction coupling intrinsic and collision degrees of freedom. The optical potential, 
which is diagonal in channel space, accounts for the average interaction between projectile 
and target. 

Usually the solution of Schrodinger's equation 

iy*(r,0 = -E*(r,0, (4) 

where E is the collision energy in the center of mass frame, is expanded as 

v|>(r,O=£^(r)0 Q (O, (5) 

a 

where ip a {v) describes the relative motion in channel a. Substituting this expansion in 
Eq. (fD we obtain the coupled channels equations (see e.g. Ref. [[HI). 

(E a - H a ) V> Q (r) = £ V a0 (r) ^(r) . (6) 



Above, E a = E - e a and H a = T + U° pt {r), where 

jjopt = yopt _ . w opt (?) 



is the optical potential in channel a. The imaginary parts have the purpose of accounting 
for the flux lost to channels neglected in the expansion of Eq.(|5]). The channel coupling 
potentials, in Eq.(0) are given by 

V Q/3 (r) = j^C(0^(r,0<MO- (8) 

A consequence of the non-Hermitian nature of H (see Eq.(|^)) is that the continuity 
equation breaks down. This can be checked following the usual procedure to derive the 
continuity equation. For each a, we evaluate ip* a {j) x [Eq.(|6])] — [Eq.(|6])]* x ijj a {r) and then 
sum the results. Assuming that V Q/ g is hermitian, we obtain 

v-£i« = l£^WIV'«(r)lVo 

a "'a 

Integrating the above equation inside a large sphere with radius larger than the interaction 
range and using the definition of the absorption cross section, we obtain the useful relation 
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^ = 4e<^i w a \^ a ) . (9) 



A. Polarization potentials 

In some coupled channel problems, it occurs that one is only interested in the elastic 
wave function. One example is the study of complete fusion in collisions involving nuclei 
far from stability, where the breakup threshold is very low. An extreme example is n Li, 
which has no bound excited state. In such cases, the coupled channel problem involves 
only the elastic and the breakup channels. Since the breakup channels contain at least three 
fragments, their contribution to complete fusion is expected to be negligible. Therefore, only 
the elastic wave function is required for the calculation of the complete fusion and breakup 
cross sections. 

In such cases, the polarization potential approach becomes very convenient. It consists 
of replacing the coupled channel equations by a single Schrodinger equation for the elastic 
state. This equation contains a polarization terra, U p ° l , added to the optical potential and its 
solution is identical to the elastic wave function obtained from the coupled channel equations. 



According to Feshbach |Tj| , the polarization potential is obtained through elimination of the 



coupled channel equations for excited states and it is given by 

U** = (MPvQd&QvPfa) . (10) 

Above, P = |0o)(0o| is the projector on the elastic channel, Q — 1 — P — J2 a ^o 1 0a) (0a I > 
and the propagator Gqq is defined as 

G ^ = E-QH Q + ze " (U) 



The wave function is then obtained by solving 

(E - H - U**)\1*) = , (12) 

which, in the position representation is written 

~E-T- U opt (r)] Vf(r) - / U pol (r, r') V(r')dV = , (13) 

where U pol (r,r') is the nonlocal potential 

U pol (v, r') = £ Vba(r) GW(£ a ; r, r') y o0 (r') . (14) 

a 

In principle, evaluating the polarization potential is nearly as hard as solving the coupled 
channel equations. However, for practical purposes it is replaced by trivially equivalent local 
potentials, which are calculated with approximations (TD],[TT|. 

B. Fusion and breakup cross sections 

With the introduction of the polarization potential, any flux going away from the elastic 
channel is treated as absorption. The sum in Eq. (|9]) is then reduced to a single term, the 
one with a = 0. The imaginary part of the potential is (henceforth we drop the superfluous 
index a, since only the elastic channel appears), 

W = W opt + W po1 , (15) 

the absorption cross section can be split as 

<?a = <?F + °bu ■ (16) 

Above, 

a F = ^Jd 3 rW opt (r) |^(r)| 2 (17) 

is identified with absorption through complete fusion and 

a bu = ^Jd 3 rW pol (r) |^(r)| 2 (18) 

corresponds to the loss of flux through the breakup channel. It includes the breakup cross 
section and also a cross section for absorption in the breakup channels, probably incomplete 
fusion. However, since for weakly bound nuclei the range of W po1 is much larger than that 
of W opt , we neglect this contribution and use the notaton a bu in Eq. (|1"8|). 
It is useful to consider the expansion in partial waves of the wavefunction, 

^E^M^)> (19) 



where k = J2fiE/h and the Ui(k,r) are solutions of the radial equation, 



h 2 Id 2 1(1 + 1) 
2J1 



dr 2 r 2 



m{k, r) + U opt (r) u,(Jfe, r) = E u t {k, r) 



normalized such that 



u,(Jfe,r -> 00) = ^ [#j (_) (^) - SiH^'(kr) 



r(+), 



(20) 



(21) 



Using the partial wave expansion in Eg. (|TT|) , the fusion cross section may be rewritten as 



a F = -E(2/+l)T/ 



where the transmission coefficient is given by 

4k <-°° 



A If POO 

Tf = l- \S h \ 2 = — dr W°P\r) \ Ul (k,r)\ 2 . 
hi Jo 

Proceeding similarly with Eq.fll^), we get 



°>» = ii £(« + i)zj 



bu 



with 



T, 



bu 



4k 

~E 



drW po \r) \ui{k,r)\ 2 . 



(22) 



(23) 



(24) 



(25) 



III. APPROXIMATIONS 



In what follows, we study different approximations to the coefficients Tf and T\ u . In order 
to fix ideas, we consider a n Li beam incident on a 12 C target, using an optical potential 
jjopt _ yopt _ j -yyopt p arame terized in the standard way: 



V opt (r) = V n (r) + V°(r 



rC, 



with the nuclear part given by 

V N (r) 

and the Coulomb one by 

V c (r) = Z p Z t e 2 /r; 

Z p Z t e 2 /2R c 



yN 
v 



1 + exp [(r — R T ) fa r ] 



for r > R, 



c 



3 -^ 



for r < Rq 



(26) 



(27) 



(28) 



Above, Z p , A p (Z t , A t ) are the atomic and mass numbers of the projectile (target), Re is the 
radius of the nuclear charge distribution, and R r is given by 



R r = r° r (A 1 J 3 + A 1/: ''\ 



< ; • 



The imaginary part is similarly parameterized as 

W° pt 



W op \i 



1 + exp [(r - Ri) /oj] ' 
with Ri defined similarly to Eq . (|29"D . We take the following parameter values: 



V° pt = -60 MeV; r° = 1.25 fm; a° r = 0.60 fm : 



(29) 



(30) 



(31) 



Wo pt = 60 MeV; 



T: 



1.00 fm; a° = 0.60fm. 



(32) 



Note that, since W opt corresponds exclusively to short range fusion absorption, r° is 
appreciably smaller than r°. 

In order to review the standard approximations in the optical potential calculations, we 
initially disconsider the breakup channels. In the absence of breakup, the imaginary part of 
the nuclear potential has a short range, and therefore fusion may be approximately described 
through an infinitely absorbing imaginary potential with a well defined radius Rf- In this 
case Tf may be estimated by TJ, the transmission coefficient through the effective potential 



Vi(r) = V opt (r) + 



h 2 1(1 + 1) 



(33) 



If one approximates the region around the maximum of VJ by a parabola, then one obtains 
the Hill- Wheeler expression for Tf pjl 



T? 



T,^T, 



HW 



1 + exp 



2n 



Bi-E 

Tiui 



-i 



(34) 



where Rb is the position of this maximum, B\ its value, and uii the curvature of Vi at r = Rb, 

1/2 

(35) 



Uui 



tf_ \d 2 Vi[r) 
ix 



dr 2 



1Kb, 



In Fig. 1, we show an example of a cross section calculated within the Hill- Wheeler 
approximation (dashed line) compared with the exact quantum mechanical calculation (full 
circles). One notices that the approximation is excellent at energies above the Coulomb 
barrier, E > Vb = Bi =0 , but worsens rapidly for E « Vb- 

The problem at low energies may be improved using the WKB approximation. The 
transmission factor is then given by 



Ti w exp (-2$) 



(36) 



fout 



where 

$ = Im \ I k(r)dr \ . (37) 

Above, 



k(r) = i yV [£ - V5001 (38) 

and r in e r ouf are the inner and outer classical turning points for the potential Vi, determined 
through the condition Vi(r in ( out j) = E . However, this approximation is not good at energies 
E ~ Bi] for E = B\ it yields 7} = 1 instead of the correct value T\ = 1/2, and even worse, 
it does not predict reflections above the barrier. Improvement is obtained by substituting 
the approximation of Eq. (|36|) by Kemble's expression |2(J below the barrier while keeping 
Hill- Wheeler's approximation above it, 

T t = (1 + e 2 *)- 1 (E < B t ); T t = T t HW (E > B t ). (39) 

We have employed this approximation in Ref. 0; it is equivalent to employing the Hill- 
Wheeler formula for all energies, albeit with the modification 

*( Bl -E) ^$; for E> B t . 

TWJi 

The cross section obtained within this approximation is depicted in Fig. 1 (solid line). We 
see that it reproduces the full quantum calculations for all collision energies. 

Let us now consider the inclusion of the breakup channels. As we have seen, this may 
be done through the introduction of an appropriate polarization potential. Such potentials 
were studied in Refs. [7U§, for pure nuclear coupling, and in PJlTj for the electromagnetic 
coupling. In |7||| only the imaginary part of the polarization potential was calculated. 
Since the real part of the polarization potential reduces the height of the potential barrier, 
this effect was simulated by a shift in the collision energy in the calculation of T\. Namely, 

T t (E) -^Ti(E + AE); AE = - V pol (R B ) ■ 

As we will see, the real part of the polarization potential plays a very important role at 
energies below the Coulomb barrier. In the case of n Li + 12 C, the breakup process is 
dominated by the nuclear coupling. Therefore we write 

wr « (r)= ^y (40) 

1 + exp [[r - Rpoi) /a] 

where R po i may be approximated by the optical potential radius, and the diffuseness a is 
given in terms of the breakup threshold energy B^ u as 

/2fibuBbu\ /,,v 



Above, fibu is the reduced mass of the fragments produced in the breakup process. In the 

?6u 



case of n Li, Br„. = 0.2 MeV and thus a = 6.6 fm. 
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The strength of the polarization potential varies with / and Ecm, and, for the partial 
waves relevant to the fusion process, is of the order of 1 MeV in the region around r rs R po i- 
Since in this work we are not concerned with its derivation, but with the approximations 
employed in the determination of the cross section, we shall adopt the constant value 

W%°\l, E) = Wl° l = 2.0 MeV . (42) 

Since the real part of the polarization plays a very important role at energies below the 
Coulomb barrier, we shall include it here. In the calculations of Andres et al. [JTT] the real 



and imaginary parts of the polarization potential have qualitatively the same strengths. For 
simplicity we then take them to be equal, i.e. 

Vg°\l, E) = Vg° l = -2.0 MeV . (43) 

The effect of the real and imaginary parts of the polarization potential are shown in 
Fig. 2. As it could be expected, the real part leads to a substantial increase in the fusion 
cross section, most evident at energies below the Coulomb barrier. On the other hand, the 
imaginary part reduces the cross section both above and below the barrier. When both the 
real and imaginary parts are included, there is a competition between the effects of the real 
and imaginary parts. With the polarization strength values considered above, suppression 
dominates above the barrier and enhancement below it. This situation was also encountered 



in the coupled channels calculations of Breitschaft et al. |L3[ and Hagino et al. [|14[ . 

The presence of a long-ranged absorption requires the introduction of modifications in the 
approximations to Tf . Now the flux that reaches the strong absorption region is attenuated 
not only because of the reflection at the barrier, but also because of its absorption into the 
breakup channel. In Ref. || it was proposed the approximation 

Tf ^T L (E + AE) ■ Pr rv , (44) 



where Ti(E+AE) is the WKB transmission factor (Eq. (pq)) evaluated at the energy E+AE 
and Pj* urv is the breakup survival probability. Within the WKB approximation we may take 



exp 



W pol (r) , 
■or 



hJ vAr) 



(45) 



where vi(r) is the local radial velocity along a classical trajectory with angular momentum 
HI. A more formal justification for Eq (pH), based on a WKB calculation with three turning 



points was presented in Ref. ||. This approximation is consistent with the results of Fig. 2. 
The enhancement due to V po1 is incuded in 7} while the suppression arising from W po1 is 
contained in p^ urv . 

In order to estimate p* urv one needs to define the classical trajectories to be employed 
in the calculation. In Ref. |J we considered pure Rutherford trajectories, neglecting the 
nuclear potential diffractive effects. These trajectories present a single turning point. The 
corresponding fusion cross section is shown in Fig. 3 as a thin line with solid circles. This 
figure also depicts the full quantum mechanical results (thick solid line). We see that al- 
though the approximation obtained with the Rutherford trajectory is reasonable at high 
energies, it breaks down at energies close and below the Coulomb barrier (Vb = 2.67 MeV). 



The inclusion of the nuclear potential in the trajectory calculations improves considerably 
the results (thin line with stars). In this case we may have, depending on the partial wave 
and collision energy, one or three turning points. The treatment with three turning points 
is not accurate in the region around the Coulomb barrier, and that is the reason why there 
are large deviations in the approximated fusion cross section. Later we will show how one 
may improve this approximation, but let us first briefly consider the breakup cross section. 
In Ref. H the breakup was calculated by considering that 



T, 



bu 



-i psurv 



(46) 



This approximation is based on the notion that Jf" corresponds to the probability of non- 
survival to the breakup process. The results depend strongly on the classical trajectory 
considered. In Fig. 4 we compare the exact quantum mechanical breakup cross section to 
the ones obtained using Eq. flJBp with different trajectories. The results are far from satis- 
factory. In particular, when the nuclear potential is included in the trajectory calculations 
the low energy breakup cross section has a completely wrong behavior. The reason for this 
discrepancy has been discussed by Takigawa et al. [§] and will be considered in further detail 
later in this section. 

Let us now develop an improved WKB approximations for Tf and Tf u . In order to 
explain them, it will be useful to rewrite the Tf coefficients in a different way. In the WKB 
approximation, the radial wave funcions with incoming (— ) and outgoing (+) boundary 
conditions are given by 



u 



(±)< 



.4 



k(r) 



exp 



±i dr k(r) 



(47) 



where 



k(r) 



n\ 



2// 



E - C/°P*(r) 



h 2 



2/ir 



-1(1 + 1)- UP° l (r) 



(48) 



The value of Tf is given by the ratio between the probability density current that reaches 
the strong absorption region, j^~\r = Rp), to the incident one j^~\r = oo), where the radial 
currents are 



J {±) (r) 



h 

2/ii 



where 



u\ (V 



du\ \r) \ (±) 
, I - u) >(r) 
dr 



From Eqs.fl47D to (PI), we obtain 



T 



3 { -\r = R F ) 

j(-)(r = oo) 



cxp 



-2$ 



.(±) 



du\ (r) 
dr 



(49) 



(50) 



$ = Im \ I dr k(r s 
Ir f 



(51) 
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If one does not include the polarization potential, the integrand in the equation above is 
real on the whole classically allowed region (note that W opt (r > Rp) = 0). In this way, only 
the classically forbidden region contributes to attenuate the current that reaches the fusion 
region (r < Rf), i-e- 

I-Tout 

$ _► $ = / drk(r), (52) 



' in 



where r in and r out are the inner and outer turning points. In this case Tf reduces to the 
expression given in Eq. Q55|). 

However, if there is long-ranged absorption as a result of the coupling to the breakup 
channels, the integrand in Eq. fl5l|) becomes complex in all the integration region. The 
contributions to the integral that defines $ from the classically allowed and forbidden regions 
may be calculated separately. In this case, Tf is written as the product of factors resulting 
from each of them. Disregarding the imaginary part of U pol (r) in the classically forbidden 
region, the corresponding factor reduces to the WKB tunneling probability TJ. On the 
other hand, in the classically allowed regions k(r) can be calculated in an approximate way. 
Assuming that the imaginary part of U pol (r) is small in comparison to the remaining terms 
in the square root appearing in Eq. (fi8|), we may take a series expansion to the lowest order, 

W pol (r) 
k {r) ^ k {r) + i — V, (53) 



%v(r) 



where 



, . . 1 
and v(r) is the local velocity, 



h 2 

E - U° pt (r) -1(1 + 1) - V p °Hr) 

2fir 2 



(54) 



v(r) = . (55) 

Since ko(r) does not attenuate the incident probability current, we obtain the same factor 
Pf urv as before. 

In our procedure we do not explicitly distinguish between classically allowed and forbid- 
den regions, and calculate $ directly from Eq. ([51]), without any of the additional approx- 
imations mentioned in the previous paragraph. In Fig. 5 we show the fusion cross section 
obtained within this approximation, compared with the exact results and with the old ap- 
proximation. We see that the present approximation yields excellent results in all energy 
regions, including the one around the Coulomb barrier where the old approximation totally 
failed. 

As noted by Takigawa et al. ||, the relationship between T^ and P^ urv that appears in 
Eq. (|46| ) is not actually correct. The reason for this is that when we calculate the survival 
probability we consider only the incident branch of the trajectory, along which the system 
approaches the strong absorption region. However, the breakup process may take place 
both on the entrance or exit branches. Let us first consider the calculation of Ref. 0, which 
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determines P^ urv along a Rutherford trajectory. The survival probability associated with 
Tf" is the one calculated along the whole trajectory, i.e. along both branches A and B in 
Fig. 6a, and not just along branch A, as it was done in the calculation of P£ urv . Since the 
contribution from both branches to the integral that defines $ (Eq. ([51])) are equal, the 
breakup probability amplitude may be written as 



ibu 



(Pi 



(56) 



If we now take into account the effect of the nuclear potential on the classical trajectory, 
the situation changes very much. For low partial waves, where E > Bi, the infinite absorp- 
tion condition in the strong absorption region allows for only an ingoing branch. On the 
other side, for partial waves for which E < B\ we may have two classical turning points, as 
illustrated in Fig. 6b. In that case all segments A, B, C do contribute to the breakup cross 
section. In this case the amplitude 7] is given by 



->6u 



,•(-) 



H/ 



jr(oo)- 3 r(r out ) 



Jl \ r in) 



;(-) 



jf'iRi 



+ 



ji +) (r out ) 



ji +) (°°; 



jf'(oo) 



(57) 



The first term in the numerator corresponds to the contribution to the breakup channel 
along incoming branch A in Fig. 6b. The second term corresponds to the other incoming 
segment, C, while the third one is the contribution associated to the exit branch B. The 
currents in this equation are given by 



3i (rout) 

Jl \ r in) 

t\R F ) 

j\ + \rout) 
jf\™) 



e" 2 * 1 j\'\oo) 



Ti Jl 



(-)/ 



out/ 



-2*2 „•(-) 
Jl 



;(-) 



Ti) j\-\r out ) 



e- 2 ^jl + \r out ), 



(58) 



where $i e $2 are given by 
$! = Im 



dr k(r) 



Tout 



$9 = Im 



Substituting the density currents in Eq. (|57|), we obtain 



dr k(r) > . 



(59) 



T, 



bu 



1-e 



-2$i 



+ e 



-2<J>i 



T, 1-e 



-2$; 



+ (1-71) (1 



,-2*] 



(60) 



The breakup cross section calculated using Eqs. (|56|) (dashed line) and (|60D (full line) are 
shown in Fig. 7, where they are compared to exact results (solid circles). We notice that the 
two approximations lead to similar results, and both are reasonably close to the exact values. 
Comparing the two curves we reach two important conclusions. One is that the inaccuracy 
in the results in Fig. 4 is due to the omission of the exit branch in the trajectories. The other 
is that in the present case nuclear effects on the trajectory are not very relevant. This is 
because the most important contributions to the breakup cross section arise from the high-/ 
partial waves. While for the energy range considered the fusion cross section converges for 
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/ = 10, the breakup one requires the inclusion of partial waves as high as I ~ 80. In this 
way, for most partial waves relevant for the breakup calculation the external turning point 
is placed outside the nuclear potential range. The situation changes somewhat when a more 
realistic potential is considered. In that case, its intensity decreases at high I values, and 
the breakup cross section becomes more sensitive to low partial waves. 

IV. CONCLUSIONS 

We have investigated the validity of commonly used approximations for complete fusion 
and breakup transmission coefficients in collisions of weakly bound projectiles at near barrier 
energies. They were tested in the concrete example of the of a n Li + 12 C collision. For the 
calculations, we adopted a typical strong absorption optical potential and, for simplicity, 
a schematic polarization potential, consistent with theoretical predictions available in the 
litterature |I§[nj. 

We have shown that the factorization of the complete fusion transmission coefficient as 
a tunneling factor times a survival probability @J§-0 may be a reasonable approximation, 
depending on the classical trajectory used in the evaluation of the latter. However, it 
is always inaccurate in the neighbourhood of the Coulomb barrier. An improved WKB 
approximation for Tf was shown to lead to very accurate values of the complete fusion cross 
section in the whole energy range of our study, both above and below the Coulomb barrier. 

We have also shown that the breakup transmission coefficient can be obtained in terms 
of survival probabilities provided that the emergent branch of the classical trajectory is 
included in the calculation. However, the accuracy of this approximation is worse than that 
for complete fusion. 
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Figure Captions 

• Figure 1: Hill- Wheeler and WKB approximations to the fusion cross section. The 
vertical arrow indicates the position of the Coulomb barrier. See text for further 
details. 

• Figure 2: Fusion cross section with different contributions of the polarization potential. 
See text for details. 

• Figure 3: Fusion cross sections obtained with different approximations employed in 
previous publications. The solid line indicates exact quantum mechanical calcula- 
tions and the remaining ones are obtained with survival probability approximation 
(Eq. Q4"3|)). The solid circles were obtained with Rutherford trajectories while the 
starts were obtained with classical trajectories taking into account both the Coulomb 
and the nuclear potentials. 

• Figure 4: Exact breakup cross section (solid line) and cross sections approximated by 
Eq. (|46f) . The solid circles were obtained with pure Rutherford trajectories while the 
stars takes into account nuclear potential effects on the trajectory. 

14 



Figure 5: Exact complete fusion cross section (solid circles) compared to the old ap- 
proximation, depicted also in Fig. 3 (stars), and with the improved WKB approxima- 
tion (solid line). See text for more details. 

Figure 6: Branches (A, incoming; B, outgoing) of the collision trajectory that con- 
tribute to the breakup process, (a) pure Rutherford, and (b) including nuclear poten- 
tial effects. In this later case, the incoming branch has an additional segment (C). 

Figure 7: Exact calculations of the breakup cross section (solid circles) compared to 
WKB calculations taking into account all branches of the classical trajectory (A, B, 
C in Fig. 6b and Eq. (|6"0D) (solid line) and taking only branches A and B in the 
Rutherford trajectory (Fig. 6a and Eq. (|56"D) (dashed line). 
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